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SECTION  I 


INTRODUCTION 

In  the  field  of  airplane  design,  there  exists  a class  of  problems  in 
which  the  motion  of  an  aeroelastic  system  may  vary  in  an  arbitrary  manner 
with  time,  depending  upon  the  Initial  conditions,  the  character  of  the  applied 
forces,  and  dynamic  response  properties  of  the  system.  Rapidly  applied 
external  forces  may  rise  from  many  sources  during  the  life  of  a modern  air- 
plane. One  of  the  most  important  disturbing  forces  is  that  produced  by  gusts. 
Detailed  study  of  this  phenomenon  In  subsonic,  transonic  and  supersonic  flows 
is  Important. 

A time-dependent  response  study  requires  simultaneous  Integration  of  the 
structural  and  aerodynamic  equations  in  time.  In  the  subsonic  and  supersonic 
cases,  the  governing  flow  equations  are  linear  and  aerodynamic  forces  depend 
upon  the  body  motion  in  a linear  fashion.  Many  unsteady  aerodynamic  theories 
are  available  fo;  the  solution  of  the  aeroelastic  problems  for  these  cases. 
Several  classical  examples,  Illustrating  the  dynamic  response  phenomena  for  the 
subsonic  case  may  be  found,  for  example,  in  Reference  1. 

On  the  other  hand,  governing  equations  for  flows  in  the  transonic  regime 
are  nonlinear  and  they  are  characterized  by  the  presence  of  shocks  on  the  air- 
foil. Because  of  these  complexities,  the  study  of  the  aeroelastic  phenomenon 
in  the  transonic  regime  Is  fairly  Involved.  From  earlier  studies  it  has  been 
observed  that  for  the  case  of  flow  over  an  airfoil  In  a free  stream  at  Mach 
numbers  near  one,  small  amplitude  motions  can  cause  large  variations  In  the 
aerodynamic  forces  and  moments.  In  addition,  phase  differences  between  the 
flow  variables  and  the  resultant  forces  may  also  be  large.  Because  of  these 
special  characteristics  of  the  transonic  flows,  the  probability  of  encountering 


aeroelastlc  Instabilities  Is  higher.  Detailed  Investigations  of  this  problem 
are  of  practical  significance. 

The  numerical  methods  for  the  computation  of  aerodynamic  forces  of  small- 
disturbance  transonic  flows  about  oscillating  airfoils  and  planar  wings  have 
been  rapidly  developed  In  recent  years.  A brief  mention  of  such  developments 
containing  32  references  can  be  found  In  Reference  2.  Based  on  these  develop- 
ments many  aeroelastlc  computations  have  been  carried  out.  A comprehensive 
state-of-the-art  review  of  the  developments  In  transonic  aeroelastlclty  Is 
given  by  Ashley  In  Reference  3. 

The  first  significant  application  of  unsteady  transonic  computations  to 
an  aeroelastlc  response  study  was  by  Ballhaus  and  Goorjlan  (Reference  4).  In 
that  study  they  Illustrated  the  use  of  the  time-integration  and  the  Indlclal 
approaches  for  the  aeroelastlc  response  problems.  They  performed  an  aeroelastlc 
response  analysis  for  a NACA  64A006  airfoil  with  a single  pitching  degree  of 
freedom.  Pitching  axis  was  assumed  to  be  at  the  mid-chord.  The  time-history 
responses  were  computed  by  using  their  computer  program  LTRAN2  for  unsteady 
transonic  flow  coupled  with  an  Integration  procedure  for  the  simple  differential 
equation  of  motion  for  the  airfoil.  The  motion  of  the  airfoil  In  a flow  at 
Mm  * 0.88  was  forced  for  the  first  few  cycles  with  an  amplitude  of  oscillation 
of  0.5  degrees  until  the  pitching  moment  became  periodic,  after  which  the  air- 
foil motion  and  the  aerodynamic  response  were  left  free  to  drive  each  other. 

As  the  structural  damping  coefficient  was  varied  parametrically  the  responses 
were  obtained  for  the  highly  damped,  neutral,  slightly  divergent,  and  highly 
divergent  cases.  It  was  shown  that  In  order  to  obtain  a flutter  solution  for 
their  example  system  at  Mm  _>  0.88,  It  Is  necessary  that  the  nonlinear  aerodynamic 
equations  be  used  and  that  the  moment  variation  lead  the  motion.  Here  It  may 
be  noted  that  linear  flat  plate  equations  can  not  be  used  to  predict  a phase 
lead  andean  not  be  used  to  obtain  a flutter  solution. 
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In  Reference  4 another  case  with  large  Initial  amplitude  of  1.5°  and 
= 0.87  was  also  considered.  This  example  illustrated  a non-sinusoidal 
pitching  moment  behavior  that  could  result  from  the  large  shock  wave  excursions 
encountered  at  larger  airfoil  motion  amplitudes. 

In  Reference  5,  Rlzzetta  performed  an  aeroelastlc  response  study  of  a NACA 
64A010  airfoil  by  simultaneously  Integrating  the  LTRAN2  aerodynamics  program 
and  the  structural  equations  of  motion.  The  study  Included  a single  degree 
of  freedom  case  with  airfoil  pitching  at  quarter  chord  and  a three  degree  of 
freedom  case  with  airfoil  pitching,  plunging,  and  aileron  pitching. 

Mach  numbers  considered  In  Reference  5 were  0.72  and  0.82,  A procedure 
similar  to  that  discussed  In  Reference  4 was  employed  to  analyze  the  single 
pitching  degree  of  freedom  case.  For  the  three  degree  of  freedom  case,  results 
were  obtained  as  plots  of  nondlmenslonal  pitching  displacement  (and  moment), 
plunging  displacement  (and  lift)  and  aileron  displacement  (and  moment)versus 
time  for  various  arbitrary  values  of  the  airfoil-air  mass  ratios.  The  other 
aeroelastlc  parameters  were  also  selected  arbitrarily.  All  the  displacement 
and  force  quantities  In  the  figures  were  divided  by  a'(0),  the  Initial  pitching 
velocity.  Dimensional  values  for  a'(0)  corresponded  to  1,3,  and  5 degrees.  By 
varying  the  values  of  y , converging  and  diverging  curves  were  obtained. 

It  was  pointed  out  In  Reference  5 that  no  attempt  was  made  to  obtain 
response  curves  corresponding  to  flutter  condition  for  the  three  degree  of 
freedom  system.  The  accuracy  of  results  was  tested  against  the  homogeneous 
solutions  obtained  for  the  structural  system  without  the  aerodynamic  forces. 

It  was  achieved  by  assuming  large  values  for  y.  At  large  values  of  y the 
quantities  In  the  aerodynamic  matrix  have  negligible  effect  on  the  response 
solution. 

Tne  aeroelastlc  response  studies  made  in  the  transonic  regime  so  far  have 
not  Included  the  single  degree  of  freedom  case  of  airfoil  plunging  and  the  two 
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degree  of  freedom  case  with  airfoil  pitching  and  plunging.  Although  the 
single  degree  of  freedom  systems  do  not  represent  the  "real  world"  problems, 
some  studies  of  them  may  give  insight  to  other  practical  aeroelastlc  problems. 

On  the  other  hand  the  earlier  studies  have  shown  that  primary  flexure- tors  Ion 
flutter  corresponding  to  the  two  degree  of  freedom  system  Is  the  most 
destructive  of  all  aeroelastlc  Instabilities.  The  "transonic  dip"  phenomenon 
reported  In  References  6,7  and  8 showed  that  flexure-torsion  flutter  may 
be  more  hazardous  at  Mach  numbers  near  one.  Studies  conducted  so  far 
(References  2, 3,6, 7, 9)  Included  only  the  results  corresponding  to  the  neutral 
stability  (flutter)  point  for  transonic  Mach  numbers.  The  response  studies 
In  the  neighborhood  of  these  transonic  neutral  stability  points  are  also  of 
significant  Interest. 

In  this  study  aeroelastlc  response  analyses  were  performed  for  a flat  plate 
and  a NACA  64A006  airfoil.  Several  cases  with  single  degree  of  freedom  of 
pitching,  single  degree  of  freedom  of  plunging,  and  two  degrees  of  freedom 
of  pitching  and  plunging  were  considered  In  detail.  Aerodynamic  forces 
required  In  this  analysis  were  obtained  by  time-integration  method  with  the 
use  of  the  LTRAN2  computer  code. 

To  evaluate  the  linear  solution  obtained  by  using  LTRAN2,  a flat  plate  was 
first  analyzed  at  ■ 0.7.  Converging  type  response  was  obtained  for  a flat 
plate  pitching  at  mid-chord  by  simultaneously  Integrating  LTRAN2  with  the 
corresponding  structural  equation  of  motion.  For  the  same  case,  response 
was  also  obtained  by  Integrating  the  aerodynamic  equation  given  by  the  quasi 
steady-state  theory  with  the  structural  equation  of  motion.  Responses  for 
pitching  angle  and  pitching  moment  were  obtained  by  using  both  methods  and 
the  two  sets  of  results  were  compared  with  each  other. 

An  example  Identical  to  that  discussed  in  Reference  4 was  also  considered. 
In  the  example  a NACA  64A006  airfoil  with  single  degree  of  freedom  pitching 


about  mid-chord  at  ° 0.88  was  considered.  In  this  study  the  aerodynamic 
coefficients  required  for  the  analysis  were  obtained  by  time  Integration 
method,  whereas  In  Reference  4 the  Indicia!  method  was  employed  to  obtain  the 
same  aerodynamic  coefficients.  Response  results  obtained  from  both  studies 
for  neutrally  stable,  diverging,  and  converging  cases  are  compared  and 
discussed.  This  example  was  also  used  to  evaluate  the  numerical  Integration 
scheme  employed  In  the  present  analysis. 

Another  case  of  a single  plunging  degree  of  freedom  system  for  a flat 
plate  at  Mach  number  of  0.7  was  also  considered.  Converging  type  plunging 
response  results  were  obtained  by  using  both  the  aerodynamic  forces  computed 
by  LTRAN2  and  quasi  steady-state  aerodynamic  theory.  The  two  sets  of  results 
were  compared. 

Finally,  the  important  case  of  a two  degree  of  freedom  system  with  plunging 
and  pitching  was  analyzed.  Two  configurations,  one  In  the  form  of  a flat 
plate  and  the  other  a NACA  64A006  airfoil  were  considered. 

The  flat  plate  was  assumed  to  pitch  about  the  mid-chord  at  Miu  - 0.7. 

Flutter  speeds  and  corresponding  reduced  frequencies  were  first  obtained  by 
varying  the  airfoil-air  mass  ratio  u for  selected  values  of  plunge  to  pitch 
frequency  ratio  w^/u^,  position  of  mass  center  X(j,  and  radius  of  gyration 
Y . The  aerodynamic  coefficients  required  for  this  analysis  were  obtained  by 
both  LTRAN2  (linear)  and  a Kernel  Function  method.  For  a selected  point  on 
the  flutter  curve  (say,  a point  corresponding  to  reduced  frequency  kc  = 0.1) 
a response  study  was  carried  out.  Neutrally  stable,  stable  and  unstable 
response  results  in  the  form  of  displacements  and  forces  were  obtained  by 
varying  the  airfoil-air  mass  ratio  u.  These  results  are  discussed.  The 
effect  of  airfoil-air  mass  ratio  on  peak  amplitudes  was  also  studied. 

A study  similar  to  that  performed  for  a flat  plate  was  also  carried  out 
for  a NACA  64A006  airfoil  pitching  about  the  quarter  chord  at  M n 0.85,  The 


plots  of  flutter  speeds  and  the  corresponding  reduced  frequencies  versus 
alrfoll-alr  mass  ratio  w were  obtained  by  using  the  aerodynamic  coefficients 
computed  by  time  Integration  method.  These  results  were  compared  with  those 
obtained  by  using  relaxation  and  Indicia!  methods  In  Reference  2.  For  a 
selected  point  on  these  flutter  curves,  a response  study  was  carried  out. 


| Neutrally  stable,  stable  and  unstable  response  results  were  obtained  In  the 

[ form  of  displacements  and  forces  by  varying  the  alrfoll-alr  mass  ratio  y. 


SECTION  II 


GOVERNING  AERODYNAMIC  AND  AEROELASTIC 
EQUATIONS  OF  MOTION 

In  this  section  the  basic  equations  employed  in  the  present  work  both  for  the 
aerodynamic  part  and  structural  part  of  the  equations  of  motion  are  discussed. 

2.1  Transonic  Flow  Equations 

Many  numerical  procedures  have  been  developed  for  computations  of  transonic 
unsteady  flow  fields  around  oscillating  two-dimensional  airfoils.  A brief 
description  of  these  developments  Is  given  In  Reference  2.  In  this  report  the 
time  Integration  method  developed  by  Ballhaus  and  Goorjlan  (Reference  8)  Is  used. 

The  simplified  basic  aerodynamic  equations,  following  the  assumptions  that 
the  flow  Is  two-dimensional,  Invlscld,  transonic  (M^  a 1),  and  that  the  velocity 
disturbances  are  small  as  compared  to  the  free  stream  velocity  U,  can  be  deduced 
from  the  general  equation  of  continuity  of  gas  dynamics  as 


kc"!*tt  * 2v:*xt 
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(1) 


where  kc  * wc/U  Is  the  reduced  frequency;  Mm  Is  the  free  stream  Mach  number;  <|> 

Is  the  disturbance  velocity  potential;  V„  * 1 - M2  - (y+l)M%w;  m Is  a function 

of  M ; and  y Is  the  ratio  of  specific  heats. 

In  deriving  the  above  equation,  the  coordinate  system  Is  fixed  with  respect 

to  the  airfoil,  and  x Is  aligned  with  the  free  stream  direction.  The  flow  Is 

defined  as  locally  subsonic  or  supersonic,  relative  to  the  fixed  coordinate 

system,  for  V > 0 or  V < 0,  respectively.  A measure  of  the  degree  of  un- 
c c 

steadiness  Is  given  by  the  reduced  frequency  kc  when  the  airfoil  Is  oscillating 
periodically  with  a frequency  w In  rad/sec. 
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As  a further  simplification  of  Equation  1,  the  frequency  of  the  transonic 
flow  can  be  assumed  as  low  so  that  k e 1 - M2  3 *-<  1.  Equation  1 may 

then  be  reduced  to 
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where  t is  the  thickness-chord  ratio  of  the  airfoil. 

Equation  2 is  suitable  for  the  time  Integration  approach.  This  approach  Is 
based  on  the  finite  difference  scheme  that  integrates  Equation  2 in  time  for 
harmonic  aerodynamic  motions  until  the  transient  states  in  the  solution  disappear 
and  the  forces  become  periodic. 

Several  numerical  procedures  are  available  to  solve  Equation  2.  Among  them 
the  procedure  developed  by  Ballhaus  and  Goorjlan  (Reference  10)  based  on  the 
alternate-direction  implicit  algorithm  has  been  proven  to  be  computationally 
efficient  and  is  being  widely  used.  This  procedure  uses  a conservative,  implicit 
finite-difference  scheme  to  tlme-accurately  integrate  the  nonlinear,  low-frequency 
transonic  small-disturbance  equation  as  defined  In  Equation  2.  A computer  code 
LTRAN2  was  developed  based  on  this  procedure.  This  code  can  be  used  to  find  the 
flow  field  solutions  for  the  airfoils  with  arbitrary  combinations  of  pitch, 
plunge  and  flap  deflections. 


2.2  Quasi  Steady-State  Aerodynamic  Equations  for  Flat  Plates 

For  subsonic  flow  conditions  exact  unsteady  aerodynamic  solutions  are 
available  for  flat  plates  with  arbitrary  pitch  and  plunge  motions.  Depending 
upon  certain  assumptions  these  solutions  may  further  be  simplified.  One  such 
solution  can  be  obtained  by  using  quasi  steady-state  assumption, 

Quasi  steady-state  theory  neglects  the  Influence  of  wake  vortices  on  the 
flow.  This  Is  valid  when  the  flat  plate  is  oscillating  with  a very  low  reduced 
frequency  (kc  5 0.0). 
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For  a flat  plate  pitching  with  an  angle  a about  an  axis  located  at  a 


distance  ba^  from  the  mid-chord  and  plunging  with  displacement  h,  the  expressions 
for  quasi  steady  aerodynamic  lifting  force,  Qh,  and  moment,  Qa  are  given  as 
(see  chapter  5 of  Reference  1), 

■ - irpb2[h  4 Ua  - ba^a]-2npUb[fl  4 Ua  4 b(l/2-a^)a]  (3) 

Qa  ■ itpbZ[bahh  - Ub(1/2  - ah)i  - b2(1/8  4 a2)a]  (4) 

4 2irpUbZ(a|1  4 1 /2 ) 4 Ua  ♦ b(l/2  - a^)a] 

where  ba^  ■ dlstanco  of  the  pitching  axis  measured  from  mid-chord  In  semi-chords 
b ■ semi-chord  length  of  the  airfoil. 

The  above  equations  are  derived  assuming  that  the  flow  Is  Incompressible 
(M^  ■ 0).  The  effect  of  Mac  \ number  is  approximated  by  multiplying  the 
aerodynamic  forces  by  the  Prandtl-Glauert  number,  i$,  which  Is  defined  as 


8 ■ rzr 

In  the  present  work  It  Is  assumed  that  the  above  equations  are  valid 
for  flat  plates  oscillating  at  - 0.7  with  reduced  frequency  kc  » 0.1. 

2.3  Aeroelastlc  Equations  of  Motion  for  Single 
and  Two  Degree  of  Freedom  Airfoil  System 

In  the  present  work  two-dimensional  airfoils  In  transonic  flow  with 
single  pitching  degree  of  freedom,  single  plunging  degree  of  freedom  and  two 
degrees  of  freedom,  plunge  and  pitch,  are  considered.  Figure  1 describes  the 
sign  conventions  and  the  variables  used  In  this  study. 

2.3.1  Pitching  Degree  of  Freedom  System 

The  equation  of  motion  for  the  case  assuming  that  airfoil  Is  pitching 
about  an  axis  located  at  a distance  a^b  from  the  mid-chord,  Is 


I a + C a + K « 
a a a 


Q_ 


(6) 


where  I = mass  moment  of  Inertia  of  the  airfoil  about  the  elastic  axis;  K = 
a a 

spring  constant  corresponding  to  pitching;  C a mechanical  damping  coefficient 


corresponding  to  pitching;  Qa  is  the  aerodynamic  moment  about  the  pitching  axis, 


defined  positive  in  positive  a direction. 


A nondimensional  form  of  Equation  3 can  be  written  as 
a"  + Aja'  + Aga  * jr 


(7) 


where  the  prime  represents  the  differentiation  with  respect  to  nondimensional 


time  wt  and  u is  the  circular  frequency  of  oscillation.  A,  “ C /I  w and  A, 

I a a L 


Ka/Ia«  are  the  nondimensional  damping  and  stiffness  parameters,  respectively. 


The  aerodynamic  moment  Qa  may  be  written  as 


Q* 


9C  cm 


(8) 


where  q-  1/2pU£,  dynamic  pressure 

c ■ full  chord  length  of  the  airfoil 
p ■ free  stream  density 


cm  ■ aerodynamic  moment  coefficient. 


The  values  for  cm  are  computed  from  LTRAN2. 


Using  Equation  8 and  the  definition  kc  ■ wc/U,  Equation  7 can  be  rewritten  as 


a"  + A^a'  + A^a 
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(9) 


This  equation  Is  Incorporated  In  LTRAN2  to  obtain  the  aeroelastlc  responses. 
The  following  derivations  are  based  on  the  assumption  that  Qa  Is  obtained 
from  quasi  steady-state  aerodynamic  theory. 

For  a single  pitching  degree  of  freedom  system,  the  aerodynamic  moment  Q 

01 

for  the  two  degree  of  freedom  case  as  given  in  Equation  4 may  be  reduced  to 
the  form 
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Qn  11  npb2(-  Ub ( 1/2  - a^)a  ~ b2(l/8  + a2)a) 

+ 2npUb2(ah  + 1/2){l)a  + b(0.6  - *h)a)  (10) 

Correcting  this  for  the  effect  of  Mach  number  and  substituting  it  into 


Equation  7 give 


djUM  + ©i a 1 + f^a  ■ 0 


01) 


where  d^  ■ 1.0  + 


VAi  + 


(0.125+aJ)p 


2(0.5-ah)h  4(a*-0.25)B 

Ttt  + unr 


fl  * A2 
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u'  * — a . airfoil-air  mass  moment  of  inertia  ratio, 
npb’ 

Equation  11  is  a second  order,  ordinary  differential  equation.  From  this 
equation  the  closed  form  solution  for  pitching  response  can  be  obtained  for 
given  initial  conditions  and  other  variables. 

2.3.2  Plunging  Degree  of  Freedom  System 

If  h Is  the  plunging  displacement  measured  positive  downwards  from  a mean 
position,  the  equation  of  motion  for  the  system  can  be  written  as 


mh  + C.  fi  + K.h  a (I 


(12) 


■V  T h ” yh 

where  m » mass  of  the  airfoil  per  unit  span;  » damping  coefficient  corresponding 

to  plunging;  spring  constant  corresponding  to  plunging;  and  » aerodynamic 
lifting  force  defined  positive  in  positive  h direction. 

The  nondimenslonal  form  of  Equation  12  is  written  as 
Qk 

V 


where  the  prime  represents  the  derivative  with  respect  to  nondimensional  time  wt; 


<5  = h/c,  nondimensional  plunging  displacement;  = C^/niw,  nondimensional  damping 
2 

parameter^  = K^/niw  , nondimensional  stiffness  parameter. 

The  aerodynamic  force  may  be  expressed  as 

Qh  - - qccj  (14) 

where  c is  the  lift,  coefficient  which  can  be  obtained  by  LTRAN2,  A negative 
sign  Is  Introduced  In  Equation  14  since  lift  force  is  defined  positive  upwards 
In  the  computer  code. 

Substituting  Equation  14  In  Equation  13  and  rearranging  give 

2c , 

«'  +8,5'  + B06  - - — -r-  (15) 


where  u » — n , airfoil-air  mass  ratio 

ffpb 

This  equation  Is  Incorporated  In  the  LTRAN2  code  to  obtain  the  aeroelastlc 
responses  for  the  plunging  case. 

The  following  derivations  are  based  on  the  assumption  that  an  explicit 
expression  for  Qh  Is  obtainable  from  the  quasi  steady-state  aerodynamic  theory. 
For  a single  plunging  degree  of  freedom  system,  the  aerodynamic  force  Qh  for 
(Equation  3)  may  be  written  as 

■ - 2 ttU pbh  - irpb^h  (16) 

Correcting  the  Qh  for  the  effect  of  Mach  number  and  substituting  it  Into 
Equation  13  give 

V" + V + V ■ 0 0?) 


where  d2  - 1 + * . «2  . B,  t . fj  ■ 


Equation  17  Is  a second  order  ordinary  differential  equation.  From  this 
equation,  closed  form  solutions  can  be  obtained  for  plunging  responses  for  given 
Initial  conditions  and  other  variables. 


I 


2.3.3  Two  Dimensional  Airfoil  with  Pitch  and  Plunge  Degrees  of  Freedom 


Figure  1 shows  an  airfoil  plunging  with  a positive  displacement  h (downward) 

and  pitching  with  a positive  angle  a (with  the  nose  up)  about  elastic  axis.  The 

elastic  axis  is  located  at  a distance  a^b  from  the  mid-chord  while  the  mass 

center  Is  located  at  a distance  x b from  the  elastic  axis.  Both  distances  are 

a 

defined  positive  when  measured  toward  the  trailing  edge. 

In  this  analysis  It  is  assumed  that  the  airfoil  Is  rigid  and  the  amplitudes 
of  oscillation  are  small.  It  is  also  assumed  that  there  is  no  coupling  In  the 
mechanical  damping.  Hence  considering  the  Inertia  forces,  damping  forces, 
elastic  forces,  and  aerodynamic  forces,  the  equations  of  motion  are 

mh  + Sol  + C^h  + K^h  » ^ (18a) 

Sh+Ia+Ca+KaKQ  (18b) 

a a a a 
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where  S ■ airfoil  static  moment  about  the  elastic  axis  and  other  quantities 
have  the  same  definitions  as  given  In  the  single  degree  of  freedom  system. 
The  nondlmenslonal  form  of  Equation  18  can  be  written  as 


C"  + xaa"  + ch5'  + 
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(19a) 
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(19b) 
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where  the  prime  represents  the  derivative  with  respect  to  nondlmenslonal  time  wt; 

5 ■ h/b  Is  nondlmenslonal  plunging  displacement;  a S/mb  Is  a coefficient  to 
be  multiplied  by  the  semi-chord  length  to  locate  the  mass  center  from  the  elastic 
axis;  wh  ■ (K^/m)  ' Is  the  uncoupled  plunging  frequency;  u Is  the  frequency  of 
oscillation  of  the  system;  r^  = (I^/mb2)*  Is  the  coefficient  to  be  multiplied 
by  the  semi-chord  to  give  the  radius  of  gyration  about  the  elastic  axis;  wn=  (K^/IJ 
Is  the  uncoupled  pitching  frequency;  ^ « (C^/mu)  is  the  nondlmenslonal  damping 
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Figure  1.  Definition  of  Parameters  for  Two  Degree  of 
Freedom  Aeroelastic  Analysis 
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parameter  corresponding  to  plunging;  t = (C^/mb  w)  is  the  nondimensi onal 
damping  parameter  corresponding  to  pitching. 

Substituting  Equations  8 and  14  into  Equation  19  and  rearranging,  the  final 
matrix  equation  for  the  response  analyses  is 


where  [M],  [C],  and  [K]  are  the  mass,  damping  and  stiffness  matrices, 
respectively,  and  they  are  defined  as 

(21a) 


[C] 

[K] 
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(21b) 


(21c) 


where  U*  = U/bwa  is  the  nondimensional  flight  speed  parameter  and  c^  and  cm 
are  the  aerodynamic  lifting  and  moment  coefficients.  These  quantities  are 
directly  obtained  from  the  computer  program  LTRAN2. 

Equation  20  is  Incorporated  into  LTRAN2  to  obtain  the  aeroelastlc 
responses  for  two  dimensional  airfoils  with  pitching  and  plunging  degrees  of 
freedom. 
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SECTION  III 


RESPONSE  SOLUTION  PROCEDURE 


In  order  to  obtain  the  aeroelastlc  response  solutions  it  Is  necessary  to 
slmulataneously  integrate  the  structural  and  aerodynamic  equations.  In  the 
present  analysis  this  integration  is  carried  out  by  a numerical  procedure. 

3.1  Direct  Integration  Method 

The  aeroelastlc  Equations  9,  15,  and  20  can  be  written  In  a general  form  as 

[M]  {u")  + [C]  (u  * } + [K]  (u)  - (p } (22) 

where  [M],  [C],  and  [K]  are  the  mass,  damping,  and  stiffness  matrices,  respectively; 
{ u}  is  the  vector  of  the  displacement  degrees  of  freedom  for  the  structure;  prime 
denotes  the  derivatives  with  respect  to  nondlmensional  time  t * wt;  { p > Is  the 
vector  of  the  aerodynamic  loads.  Equation  22  is  In  the  form  suited  for  numerical 
integration. 

Many  numerical  procedures  are  available  in  the  literature  for  the  solution  of 
the  equations  of  motion  In  a form  similar  to  Equation  22.  Discussion  of  such 
procedures  as  applied  to  dynamic  problems  of  structures  can  be  found  In  references, 
such  as  11.  In  the  present  analysis  the  common  direct  integration  procedure  is 
employed  to  find  the  time-history  dynamic  responses  of  the  aeroelastlc  system. 

The  solution  Is  obtained  by  using  a step-by-step  time  Integration  finite  difference 
approach. 

Assuming  a linear  variation  of  acceleration,  the  velocities  and  displacements 
at  the  end  of  a small  time  Interval  At"  can  be  expressed  as 


(u'l-  » (u'It  - + ^fu">r  - + ™(u"l- 
1U  't  1 'l-At  t-At  2lu  't 


(23a) 
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Substituting  Equations  23a  and  23b  Into  22  yields 


where 


{u"}£  = [F]  [{p}-t  - [C]{ v)  - [K]{w>] 

(24) 

[F]  - [[M]  + 4r  [C]  + CK]]'1 

(25a) 

(V)  ■ + Y 

(25b) 

(w)  . ♦ 4- 

(25c) 

For  small  amplitude  (structurally  linear)  problems, 

matrix  [F]  need  only  be 

formed  once  since  It  Is  Independent  of  time.  On  the  other  hand  the  aerodynamic 
load  vector  { p > depends  upon  the  displacement  degrees  of  freedom  (u),  the  time 
derivatives  of  {ul,  the  reduced  frequency  k and  the  free  stream  Mach  number  M . 
The  vector  {p}  Is  obtained  by  numerically  solving  the  transonic  aerodynamic 
Equation  2 with  the  use  of  LTRAN2.  The  values  of  {u)  and  the  time  derivatives  of 
(u)  used  for  computing  { p}  are  based  on  the  values  obtained  at  the  time  t-At. 


3.2  Solution  for  Aerodynamic  Forces  {p}  by  LTRAN2 

The  basic  LTRAN2  code  employs  a non-iterative  alternating  direction  Implicit 
(ADI)  scheme  to  advance  the  solution  for  the  perturbation  potential,  4>,  from  one 
time  Interval  to  tne  next  at  each  grid  point  In  the  computational  flow  field. 
Differencing  In  x-dlrectlon  Is  of  mixed  type  which  has  been  quite  successful  In 
maintaining  stability  for  both  subsonic  and  supersonic  flow  regions.  The 
conservative  form  of  the  equation  is  preserved,  which  Is  essential  for  the 
proper  description  of  the  shock  wave  motions.  While  the  ADI  scheme  has  no  time 
step  limitation  for  stability  based  on  classical  linear  stability  analysis, 
Instabilities  may  be  generated  by  t.he  motion  of  shock  waves  due  to  the  mixed 
differencing.  Hence,  the  time  interval  a?  must  be  chosen  such  that  shock  waves 
do  not  travel  more  than  one  mesh  point  In  the  x-direction  over  a single  time 
step.  Based  on  this  procedure  pressure  distributions  at  any  time  can  be  computed. 
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By  Integrating  these  pressure  curves,  aerodynamic  lifting  forces  and 
moments  are  computed. 

At  every  time  step  LTRAN2  requires  the  effective  Induced  angle  of  attack 
and  Its  time  derivatives  as  Input.  In  general,  for  the  two  degree  of  freedom 
system  described  In  Figure  1 the  Induced  angle  of  attack  can  be  defined  as 

a1  " a + O'  (26a) 

or 

k 

« a + e'  (26b) 

a and  £'  at  every  time  step  can  be  computed  from  the  direct  Integration 
scheme  discussed  earlier  In  this  section. 

3.3  Summary  of  Step  by  Step  Time  Integration  of  Aeroelastlc  Equations 

The  procedure  discussed  In  this  section  was  Incorporated  in  LTRAN2  In  the 
following  manner.  For  a set  of  starting  values  of  (u),  {u ' > , and  {p},  (say, 
known  at  time  t-A$)  the  acceleration  vector  {u" } at  time  l was  computed  from 
Equation  24.  Based  on  the  known  acceleration  vector  {u"l,  the  displacement 
vector  {u}  and  velocity  vector  {u ’ > at  time  l were  computed  from  Equations  23a 
and  23b,  respectively.  From  these  quantities  the  effective  Induced  angle  of 
attack  a^,  and  Its  time  derivatives  were  computed  for  time  t.  This  Induced 
angle  of  attack  and  other  required  quantities  were  then  read  Into  the  LTRAN2 
code  and  the  new  aerodynamic  load  vector  {p > at  time  t was  computed.  At  this 
stage  all  the  quantities, namely , tu } , { u ' ) , { u " > , and  {p}  at  time  l were  known 
so  that  further  computations  for  time  l + At  can  be  carried  out.  This  process 
was  repeated  for  every  time  step. 
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SECTION  IV 


RESULTS 

Aeroelastlc  response  analyses  were  carried  out  for  two  airfoil  configurations, 
a flat  plate  and  a NACA  64A0G6  airfoil.  In  both  cases,  the  airfoils  were  con- 
sidered as  single  and  two  degree  of  freedom  systems. 

4,1  Finite  Difference  Computational  Grid 

The  size  and  the  pattern  of  the  grid  play  an  Important  role  In  obtaining 
accurate  results  from  LTRAN2.  Because  of  the  various  reasons  discussed  In 
Reference  10,  a fairly  fine  mesh  Is  required  in  order  to  obtain  acceptable 
solutions.  However,  the  mesh  size  Is  limited  by  the  capacity  of  the  core 
memory  avallalbe  In  a computer. 

In  this  analysis  a 79  x 99  finite  difference  mesh,  with  79  grid  points  in 
the  vertical  (y)  direction  and  99  grid  points  in  horizontal  (x)  direction,  was 
employed  for  final  computations.  Details  of  this  grid  were  kindly  supplied  by 
Ba 11 haus  and  Goorjian.  This  was  the  maximum  grid  size  that  could  be  practically 
used  on  the  CDC  6500  computer  available  at  Purdue  University.  A schematic 
diagram  of  the  portion  of  the  grid  near  the  airfoil  Is  shown  In  Figure  2. 

This  grid  has  a smooth  non-uniform  pattern  and  It  Is  symmetric  about  the 
y - 0 line.  The  spaclngs  of  the  grid  points  are  smaller  near  the  leading  and 
trailing  edges  In  x-dlrectlon  and  near  the  mean  chord  line  In  the  y-directlon. 

The  spaclngs  are  gradually  enlarged  as  the  grid  points  are  farther  away  from  the 
airfoil.  Thus  the  grid  boundaries  are  located  sufficiently  away  from  the  airfoil 
both  In  x-and  y-directlons.  From  leading  edge  to  trailing  edge,  a total  of  33 
grid  points  Is  used  so  that  an  accurate  representation  of  the  pressure  distribution 
can  be  obtained,  particularly  when  there  Is  a shock.  Other  salient  features  for 
this  grid  are  given  as  follows: 
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(1)  Minimum  spacing  In  x-directlon  = 0.0033c 

(2)  Minimum  spacing  in  y-directlon  = 0.02c 

(3)  Distance  between  the  upstream  grid  boundary  and  the 
leading  edge  = 1033.53c 

(4)  Distance  between  the  downstream  grid  boundary  and  the 
leading  edge  ■ 855.91c 

(5)  Distance  between  grid  boundary  and  the 
mean  chord  line  = 811.12c. 

4.2  Steady-State  Computations 

The  steady-state  solution  Is  required  as  an  Initial  condition  for  the  unsteady 
computations.  In  this  analysis  the  steady-state  solution  was  obtained  by  using 
the  successive  line  over-relaxation  method  (SLOR)  available  In  the  LTRAN2  code. 

In  the  process  of  obtaining  the  steady-state  solution,  multiple  grid 
computations  were  carried  out  In  order  to  accelerate  the  rate  of  convergence. 

First,  a converged  steady-state  solution  was  obtained  for  a coarse  grid  (35x38). 
This  solution  was  then  Interpolated  to  form  a starting  solution  for  a medium 
grid  (53x79).  A converged  solution  was  then  obtained  for  the  medium  grid.  This 
solution  was  again  Interpolated  to  form  a starting  solution  for  the  final  com- 
putational grid  (79  x 99). 

The  SLOR  computations  were  carried  out  for  the  fine  grid  (79x99)  for  several 
hundred  Iterations.  This  Iterative  process  was  stopped  when  the  variation  In  the 
perturbation  velocity  potential  at  all  grid  points  between  consecutive  Iterations 
reached  a value  of  about  4x10  . This  required  about  600  Iterations.  The  steady 

solution  obtained  at  this  stage  was  used  to  plot  the  pressure  curves  and  also 
used  as  a starting  solution  for  the  'ady  computations. 

4.3  Unsteady  Computations  (Forced  Motion) 

Steady-state  results  obtained  by  the  SLOR  method  were  used  as  the  Initial 
conditions  for  unsteady  computations.  In  order  to  obtain  unsteady  results  the 
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airfoil  was  first  subjected  to  a harmonic  forced  motion.  In  most  of  the  cases 
the  Induced  angle  of  attack  was  varied  sinusoidally  with  an  amplitude  of  0.01 
radians (0.574°)  and  Equation  2 was  Integrated  In  time  by  LTRAN2,  During  the 
process  of  time  Integration,  120  time  steps  per  cycle  were  used.  After  some 
duration  of  time  the  effect  of  the  Initial  conditions  on  the  unsteady  solution 
became  negligible  such  that  the  aerodynamic  force  coefficients  cfc  and  c^ 
became  periodic.  In  general,  four  to  six  cycles  of  forced  motion  were  required 
to  obtain  fairly  periodic  results  from  URAN2,  However,  the  exact  time  at  which 
the  forced  motion  had  to  be  stopped,  depended  upon  the  type  of  Initial  conditions 
specified  for  free  motion.  For  all  the  cases  studied,  forced  motion  was  stopped 
such  that  aj,  and  a'j  were  0.0,  0.01,  and  0,0,  respectively,  with 

In  the  unit  of  radians.  From  this  stage  onwards  free  motion  conditions  were 
simulated  by  simultaneously  Integrating  the  structural  and  aerodynamic  equations 
and  allowing  the  n.rfoll  motion  and  aerodynamic  response  to  drive  each  other. 

Due  to  Inherent  nonlinearity  In  the  transonic  aerodynamic  equation,  the 
free  motion  response  solution  depends  upon  the  Initial  values  for  c^,  aj,  and 
a'j.  Physically  these  Initial  conditions  may  be  used  to  represent  the  Impulsive 
forces  (gusts)  striking  the  airfoil. 

4.4  Response  Analysis  of  Single  Degree  of  Freedom  Systems 

4.4.1  NACA  64A006  Airfoil  Pitching  at  M « 0.88 

A case  of  a NACA  64A006  airfoil  pitching  about  the  mid-chord  at  ■ 0.88 
and  kc  ■ 0.1  was  considered,  This  case  was  selected  In  order  to  verify  the 
present  formulations  and  also  compare  the  results  with  those  already  obtained  by 
Ballhaus  and  Goorjlan  (Reference  4). 

The  airfoil  configuration  data  required  for  this  case  were  obtained  from 
Reference  12.  The  airfoil  configuration  Is  shown  in  Figure  3. 
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MCA  64 A 006  Airfoil  Configuration 


The  steady-state  solution  was  obtained  by  the  SLOR  method.  A plot  of  the 
steady-state  pressure  curve  Is  shown  in  Figure  4.  In  this  figure  it  may  be 
observed  that  there  is  a fairly  strong  shock  on  the  airfoil  between  x/c  ■ 0.64 
and  0.69. 

The  governing  differential  equation  of  motion  for  this  case  Is 

*"  + A^'  + Aza"  - A3cm  (27) 

p 

where  Aj  ■ u/8u'kc  (see  Equation  9). 

To  obtain  the  present  response  solution,  It  was  first  necessary  to  obtain  a 
solution  for  the  neutrally  stable  case.  The  moment  coefficient  |cm  | and  phase 
angle  (phase  lag  between  cffl  and  a)  were  obtained  by  the  time  Integration  method 
by  considering  five  cycles  of  forced  motion.  The  pitching  motion  was  sinusoidally 
forced  with  an  amplitude  of  0.5°.  The  values  of  |c  | and  t|>  obtained  from  this 
analysis  were  0.8S65  and  -67.5  degrees,  respectively.  On  the  other  hand  the 
Indlclal  method  was  employed  in  Reference  4 to  obtain  these  results.  The  cor- 
responding values  of  |cm  | and  obtained  In  Reference  4 were  0,8617  and  *68.87 
degrees,  respectively.  The  differences  between  the  present  results  and  those 
obtained  In  Reference  4 may  be  mainly  due  to  the  two  different  methods  employed. 

J a. 

By  assuming  harmonic  motion  with  frequency  w,  l.e.,  a * a(e1u  , and  by 
relating  the  aerodynamic  and  structural  constants,  Equation  27  yields 


-A3|cJ,1nt 

(28a) 

1 * A3|c,m|cos* 

(28b) 

Assuming  the  value  of  A^  equal  to  1.333  (same  as  that  used  in  Reference  4) 
and  substituting  the  values  described  earlier  for  |c  | and  g>  Into  Equations  28, 
A-|  and  A2  were  obtained,  respectively,  as  1.05  and  1.437.  These  values  correspond 
to  the  neutrally  stable  case.  It  is  noted  that  In  Reference  4,  A^  and  A2  were 
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Figure  4.  Distribution  of  Steady  Pressure  Coefficients  for  NACA  64A006 
Airfoil  at  M » 0.88  by  L TRAN 2 . 


25 


-■  gtfjhgprtgi  M heewl  km 


obtained  as  1,072  and  1.414,  respectively.  The  aeroelastlc  response  analysis  was 
then  carried  out  by  simultaneously  Integrating  the  aerodynamic  and  the  airfoil 
motion  equations.  Different  types  of  responses  were  obtained  by  parametrically 
varying  the  structural  damping  parameter  A^ . The  motion  was  forced  for  the  first 
five  cycles  after  which  the  pitching  moment  became  periodic  and  the  airfoil 
motion  and  aerodynamic  response  were  left  free  to  drive  each  other.  The  amplitude 
of  the  pitching  angle  was  0.5  degrees.  The  Initial  conditions  for  free  motion  were 
corresponding  to  a(0)  * 0.0,  and  c*'(0)  « 0.5. 

Response  was  obtained  for  values  of  the  parameters  corresponding  to  the 
neutrally  stable  case,1.e.,A^  » 1.05,  Ag  * 1.437  and  A^  ■ 1.333.  Figure  5 
shows  the  response  curves  for  the  pitching  angle  a and  the  corresponding  pitching 
moment  cn).  In  both  curves,  the  first  cycle  Is  due  to  forced  motion  and  the 
remaining  cycles  are  due  to  free  motion.  Because  the  parameters  corresponding  to 
the  flutter  solution  were  chosen,  the  response  In  free  motion  is,  as  expected, 
very  nearly  periodic. 

A diverging  response  was  obtained  when  A^  was  assumed  as  0.5.  The  results 
for  the  responses  of  both  the  pitching  angle  and  the  pitching  moment  are  shown  In 
Figure  6.  A converging  response  was  obtained  when  A^  was  assumed  as  1.5  and 
the  results  are  shown  In  Figure  7.  These  results  are  similar  to  those  presented 
In  Reference  4. 

The  above  studies  Indicate  that  the  general  formulation  of  the  problem  Is 
correct  and  the  numerical  Integration  procedure  employed  is  quite  accurate, 

4.4.2  Flat  Plate  Pitching  at  f1ro  ■ 0,70 

In  the  computer  code  LTRAN2,  there  Is  a capability  to  solve  the  linear  form 
of  Equation  2,  This  can  be  done  by  setting  the  Input  value  for  y equal  to  -1.0. 

A response  analysis  for  a flat  plate  pitching  about  mid-chord  at  Mach  number  M ■*  0. 
was  performed  for  a reduced  frequency  k * 0.1.  The  present  LTRAN2  (linear) 
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Pitching  Motion  Response  for  NACA  64A9C6  Airfoil  Pitching  about  Mid-Chord  (Neutrally 
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Figure  6.  Pitching  Motion  Response  for  NACA  64A006  Airfoil  Pitching  about  Mid-Chord 
(Diverging  Case:  A.  = 0.5,  A_  = 1.437,  A_  = 1.333). 


results  were  compared  with  that  obtained  by  quasi  steady-state  aerodynamic  theory 
which  is  valid  for  reduced  frequencies  of  the  order  kc  = 0.1. 

The  flat  plate  configuration  was  simulated  in  LTRAN2  by  using  data  corresponding 
to  a flat  plate  of  a thickriess-to-chord  ratio  of  0.01.  In  order  to  avoid  numerical 
difficulty,  the  data  for  the  airfoil  configuration  were  started  at  x/c  ■ 0.003 
rather  than  at  x/c  = 0.0  (leading  edge).  The  unsteady  computations  were  started 
with  zero  steady-state  pressure  conditions. 

Equation  9 as  given  In  section  II  is  the  governing  differential  equation  of 
motion  for  this  case.  The  aerodynamic  equation  was  integrated  in  time  for  four 
cycles  by  forcing  a sinusoidal  variation  of  pitching  angle  with  amplitude  of  0.01 
radians.  The  free  motion  was  started  at  the  end  of  the  fourth  cycle.  The  initial 
conditions  obtained  for  free  motion  were  corresponding  to  a(0)  n 0.0  and 
ct'(0)  * 0.01.  The  structural  parameters  for  free  motion  were  so  selected  that 
a converging  type  response  could  be  obtained.  The  values  for  the  damping  parameter 
A-| , the  stiffness  parameter  A2  and  the  airfoil -air  moment  of  inertia  ratio  y ‘ 
were  0.5, 1.5  and  1000,  respectively.  The  value  assumed  for  y'  was  quite  high 
when  compared  to  the  actual  values  for  aircraft  wings.  However,  this  number  was 
required  in  order  to  obtain  a response  solution  that  can  be  compared  with 
quasi  steady-state  theory. 

The  converging  type  response  curves  obtained  for  pitching  angle  a and  pitching 
moment  cm  are  shown  in  Figure  8.  In  the  same  figure  the  responses  obtained  by 
employing  the  quasi  steady-state  aerodynamic  theory  are  also  shown.  These  re- 
sponses were  obtained  by  solving  the  differential  Equation  11  for  the  same  values 
of  Aj,  A2,  y',  a (0)  and  a'(0)  as  used  for  LTRAN2.  The  two  sets  of  solutions  are,  in 
general,  in  fairly  good  agreement. 

Small  differences  in  amplitude  and  phase  angle  in  Figure  8 between  the  two 
sets  of  results  are  mainly  due  to  the  difference  between  the  two  methods.  It  may 
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Figure  8.  Comparison  of  the  Pitching  Responses  Obtained  by  LTRAN2  and  Quasi  Steady 
State  Theory  for  a Flat  Plate  Pitching  about  Mid-Chord. 
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be  observed  that  there  Is  no  phase  lag  between  the  pitching  angle  and  the  cor- 
responding moment  In  the  quasi  steady-state  solution  but  there  Is  some  phase 
lag  In  the  LTRAN2  s '•utlon.  Therefore,  the  effective  aerodynamic  dampings  In 
the  system  between  the  two  methods  are  different  which  may  cause  the  discrepancy 
between  the  two  sets  of  curves.  Such  discrepancy  becomes  larger  for  smaller 
values  of  u'  since  aerodynamic  damping  Is  Inversely  proportional  to  u1. 

At  this  point,  It  may  be  noted  that  LTRAN2  Is  based  on  the  low  frequency  ap- 
proximation which  Is  more  valid  for  Mach  numbers  near  unity.  The  discrepancy 
found  In  the  two  sets  of  curves  may  also  be  partly  due  to  the  low  value  of  Mach 
number  (0.7)  considered. 

4.4.3  Flat  Plate  Plunging  at  M ■ 0.7 

In  this  case,  response  results  were  obtained  for  a flat  plate  with  only  a 
single  plunging  degree  of  freedom.  The  Mach  number  and  reduced  frequency  k 

V 

considered  were  equal  to  0.7  and  0.1,  respectively.  Flat  plate  conditions  were 
simulated  In  the  same  manner  as  discussed  earlier.  Results  were  obtained  both 
from  LTRAN2  (linear)  and  quasi  steady-state  theory. 

The  governing  differential  equation  of  motion  for  this  case  Is  given  In 
Section  II  as  Equation  15.  The  aerodynamic  equation  was  first  Integrated  by 
LTRAN2  for  four  eyries  by  forcing  a sinusoidal  plunging  motion  with  amplitude 
of  plunging  displacement  6 ■ h/c  ■ 0.1.  This  corresponds  to  an  amplitude  of 
0.01  radians  for  the  Induced  angle  of  attack  ■ kc«! 

The  free  motion  was  started  with  Initial  conditions  corresponding  to 
a.| (0)  ■ 0.0  (6  ■ 0.1  ) and  a^(0)  * 0.01  (s'  ■ 0.0).  The  structural  parameters 
were  selected  so  that  a converging  type  response  could  be  obtained.  The  values 
assumed  for  the  structural  damping  parameter  , the  structural  stiffness 
parameter  02  and  the  alrfoll-alr  mass  ratio  y were  equal  to  0.0,  1.0,  and  100, 
respectively.  Response  curves  obtained  from  LTRAN2  for  plunging  displacement  fi 
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and  lifting  force  coefficient  c.  are  shown  In  Figure  9.  In  the  same  figure  the 
corresponding  results  obtained  by  quasi  steady-state  theory  are  also  shown. 

These  results  were  obtained  by  solving  the  differential  equation  17  based  on  the 
same  Initial  conditions  and  structural  parameters  as  used  In  L.TRAN2, 

Response  results  obtained  by  both  methods  compare  fairly  well.  The  small 
differences  In  phase  angles  and  amplitudes  may  be  mainly  due  to  the  difference 
in  the  methods. 

4.5  Two  Degree  of  Freedom  Systems  (Pitching  and  Plunging) 

In  the  response  studies  of  the  two  degree  of  freedom  systems,  two  coses  were 
analyzed.  One  case  Is  a flat  plate  pitching  about  the  mid-chord  at  M - 0.7 

oo 

and  the  other  Is  a NACA  64A006  airfoil  pitching  about  the  quarter  chord  at 
■ 0.85.  Response  results  were  obtained  by  parametrically  varying  the  airfoil- 
air  mass  ratio  p.  In  all  the  examples,  the  values  of  reduced  frequency  k„  and 
radius  of  gyration  r^  were  assumed  as  0.1  and  0.5, respectively,  The  mechanical 
damping  was  assumed  as  zero  " ca  “ 0)> 

4.5.1  Flat  Plate  Plunging  and  Pitching  about  Mid-Chord  at  M ■ 0.7 

In  this  case  response  studies  were  performed  for  a flat  plate  pitching  about 
mid-chord  at  a subsonic  Mach  number  M -0.7.  Both  the  LTRAN2  (linear)  code  and 
the  Kernel  Function  method  were  used  so  that  the  two  sets  of  results  can  be 
compared. 

Before  starting  the  response  analysis  It  was  necessary  to  find  the  aero- 
elastic  parameters  corresponding  to  a neutrally  stable  condition.  Equations  20 
discussed  In  Section  II  are  the  governing  equations  of  motion  for  this  case.  The 

it t 

aeroelastlc  parameters  Include:  flight  speed  parameter  U , airfoil-air  mass 
ratio  u,  position  of  the  mass  center  xa,  plunge  to  pitch  frequency  ratio 
radius  of  gyration  ra,  damping  parameters  and  and  position  of  elastic  axis  a^. 
The  values  for  x , ra,  r.a,  and  ah  were  assumed  as  0.0,  0.2,  0.5, 
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Figure  9.  Comparison  of  the  Plunging  Responses  Obtained  by  LTRAN2  and  Quasi 
Steady-State  Theory  for  a Flat  Plate. 
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In  order  to  obtain  the  values  of  k , U and  y correspondinq  to  the  neutrally 
stable  condition,  the  U-g  method  was  used.  A detailed  discussion  of  this 
method  is  given  In  Reference  2. 

To  use  the  U-g  method, it  Is  necessary  to  obtain  a table  of  aerodynamic 
coefficients  similar  to  those  obtained  In  Reference  2,  Table  1 shows  the  aero- 
dynamic coefficients  of  c,  , c , c and  cm  obtained  by  using  the  time  Integration 

K'd  ma 

method  of  LTRAN2  (linear)  and  Kernel  Function  method  for  various  values  of 
reduced  frequency  k_.  Agreement  Is  good  between  the  two  methods. 

A flutter  analysis  was  carried  out  for  various  values  of  alrfoll-alr  mass 
ratio  y by  the  U-g  method.  Figure  10  shows  the  curves  for  flutter  speed  and  correspond- 
ing reduced  frequency  versus  the  airfoil -air  mass  ratio  y obtained  by  both  methods. 
Excellent  agreement  Is  observed. 

Any  point  on  the  curves  shown  In  Figure  10  represents  a neutrally  stable 
state  whereas  the  zonos  below  and  above  the  curves  represent  stable  and  unstable 
conditions,  respectively.  Response  analyses  were  carried  out  for  stable,  neutrally 
stable  and  unstable  conditions  by  selecting  the  aeroelastlc  parameters  based  on 
any  selected  point  In  the  flutter  speed  curves.  Because  the  time  Integration 
mathod  was  used  In  performing  the  response  studies,  It  was  essential  that  the 
flutter  curves  be  obtained  by  ur.lng  the  time  integration  method  instead  of  the  Indicia! 
and  relaxation  methods  as  used  In  Reference  2. 

The  governing  equations  of  motion  for  this  case  are  given  by  Equation  20  of 
Section  II.  To  begin  with,  the  aerodynamic  equation  was  Integrated  In  time  for 
four  cycles  by  forcing  a sinusoidal  pitching  motion  with  amplitude  of  0.01 
radians.  In  the  fourth  cycle  the  aerodynamic  force  responses  became  almost 
periodic.  After  that,  tha  free  motion  was  started  by  simultaneously  Integrating 
the  structural  and  aerodynamic  equations.  This  was  started  with  Initial  conditions 
with  C(0)  » 0,  *'(0)  - 0,  a(0)  ■ 0,  and  a'(0)  - 0.01. 
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Table  1 


AERODYNAMIC  COEFFICIENTS  FOR  FLAT  PLATE 
PITCHING  ABOUT  MID-CHORD  AT  M - 0.7 


Aerodynamic 

loeff. 

Method 

— 

Reduced  Frequency  kc  - 

uc/U 

0. 

05 

0.10 

0. 

16 

Real 

Itnag. 

Real 

Imag. 

Real 

Imag. 

c 

1 

0.0669 

0.4225 

0.1700 

0.8000 

0.3038 

1.134 

2 

0.0616 

0.3974 

0.1666 

0.7190 

0.2719 

0.9876 

c 

9 

1 

8.449 

-1.338 

8.001 

t 

O 

7.557 

-2.025 

'a 

2 

7,964 

-1.113 

7.240 

-1.488 

6.671 

-1.570 

c 

in „ 

1 

0.0151 

0.0952 

0.0479 

0.1787 

0.0954 

0.2486 

6 

2 

0.0174 

0.0990 

0.0501 

0.1780 

0.0865 

0.2425 

cm 

m 

1 

1.904 

- .3016 

1.787 

- .4788 

1.657 

- .6361 

a 

2 

1.985 

- .3500 

1.794 

- .5117 

1.640 

- .5996 

Method  1 ! Time  Integration  by  LTRAN2  (linear). 
Method  2:  Kernel  Function  Method. 
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Figure  10.  Effect  of  Airfoil-Air  Mess  Ratio  on  Flutter  Speed  for  Flat  Plate 
at  M * 0.7 
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Aeroelastlc  parameters  for  the  neutrally  stable  condition  were  selected 
from  a point  on  the  flutter  speed  curves  In  Figure  10  at  k • 0.1.  The  cor- 
responding  values  of  u and  U Mart  equal  to  289.5  and  7.23,  respectively.  These 
values  were  substituted  Into  Equation  20  and  the  response  analysis  was  carried 
out.  Response  curves  for  pitching  angle  a end  the  corresponding  pitching 
moment  cw  for  about  six  cycles  are  shown  In  Figure  11.  After  small  initial  dis- 
turbances (due  to  Initial  conditions)  both  response  curves  show  a perfect  periodic 
behavior.  Similar  curves  obtained  for  plunging  displacement  t and  corresponding 
lifting  force  c&  are  shown  In  Figure  12.  These  curves  also  show  a perfect 
periodic  behavior  after  some  small  initial  disturbances.  It  1$  seen  that  the 
neutrally  stable  conditions  obtained  by  the  present  response  method  are  In  good 
agreement  with  those  obtained  by  the  U-g  method  based  on  both  the  time  integration 
and  the  Kernel  Function  method. 

In  Figure  11,  the  response  cut  i for  pitching  displacement  a end  pitching 
moment  dm  for  the  stable  and  unstable  conditions  ere  also  shown.  Thase  were 
obtained  by  changing  the  alrfoll-alr  mass  ratio  w.  For  unstable  response  the  » 
value  assumed  was  10X  lower  than  that  corresponding  to  the  neutrally  stable 
condition.  This  point  Is  In  the  unstable  zone  of  Figure  10.  On  the  other  hand, 
for  stable  response  the  u value  selected  was  10X  higher  than  that  corresponding 
to  the  neutrally  stable  condition.  Tho  corresponding  stable  and  unstable  responses 
for  plunging  displacement  e and  lifting  force  cfc  are  also  shown  In  Figure  12. 

It  Is  seen  that  the  stable  and  unstable  conditions  obtained  by  the  U-g  method 
produced  converging  and  diverging  responses,  respectively. 

It  was  also  of  Interest  to  study  the  effect  of  alrfoll-alr  mass  ratio 
(altitude)  w on  the  peak  amplitudes  of  the  response  curves.  In  this  analysis, 
the  peak  amplitudes  corresponding  to  the  second  free  cycle  were  considered. 

Figure  13  shows  plots  of  the  ratio  of  the  second  cycle  peak  amplitude  to  the 
amplitude  of  the  neutrally  stable  curve  versus  the  ratio  of  alrfoll-alr  mass 


wh/cjot-0.2toh- 0.0,  ,Xar  0.0 0.5 , U/bw^ 7.2327,  kc* 0.1 


Figure  11.  Effect  of  Airfoil-Air  Hass  Ratio  (Altitude)  on  Pitching  Rotation  (a)  and 
Pitching  Hcaent  (c  ) for  Flat  Plate  at  H = 0.7. 
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Figure  13.  Effect  of  Airfoil-Air  Mass  Ratio  on  Peak  Amplitudes  of 
Responses  for  Flat  Plate  at  M ■ 0.7. 
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ratio  u to  the  airfoil-air  mass  ratio  at  the  neutrally  stable  condition.  The 
plots  include  four  curves  for  plunging  displacement,  lifting  force,  pitching 
displacement  and  pitching  moment,  respectively. 

From  Figure  13  it  may  be  observed  that  the  curve  corresponding  to  plunging 
displacement  Is  steeper  than  the  other  three  curves.  All  the  curves  assume  a 
similar  trend  below  the  neutrally  stable  point.  The  curves  for  lifting  force  and 
moment  are  very  close  to  each  other  and  they  virtually  coincide.  Such  curves 
may  be  used  to  study  the  aeroelastic  behavior  of  the  system  near  the  neutrally 
stable  point. 

4.5.2  NACA  64A006  Airfoil  Plunging  and  Pitching  about  Quarter  Chord  at  M u 0.85 

In  the  flutter  analysis  performed  in  References  6 and  7,  a "transonic  dip" 
phenomenon  was  observed  for  the  NACA  64A006  airfoil.  Flutter  speeds  reached  a 
minimum  in  the  neighborhood  of  M ■ 0.85.  From  the  aeroelastic  point  of  view, 
this  Mach  number  Is  the  critical  value.  Hence,  in  this  study  aeroelastic  response 
analysis  was  carried  out  at  M ■ 0.85. 

In  this  case  a NACA  64A006  airfoil  pitching  about  the  quarter  chord  (ah“-,5) 
at  M ■ 0.85  was  considered.  The  procedure  for  obtaining  the  responses  is  the  same 
as  that  described  for  the  case  of  the  flat  plate. 

The  steady- state  pressure  distribution  curve  obtained  from  LTRAN2  (non- 
linear) is  shown  in  Figure  14.  In  the  same  figure  pressure  curves  obtained  by 
STRANS2  (Reference  2)  and  by  experimental  method  (Reference  13)  are  also  shown  for 
comparison.  In  general,  all  three  curves  agree  fairly  well. 

Using  the  steady-state  initial  conditions,  unsteady  computations  were  carried 
out  by  using  LTRAN2  (nonlinear).  The  time  integration  method  was  employed  to 

obtain  the  aerodynamic  coefficients  c^  , c^  , cm  and  cm  . The  coefficients  for 

5 a 6 a 

various  values  of  reduced  frequency  kc  are  shown  in  Table  2.  In  the  same  table,  the 
aerodynamic  coefficients  computed  by  UTRANS2  and  the  Indlcial  method  in  Reference  2 
are  also  given.  In  general,  all  three  methods  agree  well.  It  may  be  noted  that 
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the  results  obtained  by  the  time  Integration  method  are  expected  to  be  more  accurate 
than  those  obtained  by  the  other  two  methods  for  kc  < 0.2. 

Based  on  the  coefficients  given  In  Table  2,  flutter  analysis  was  carried  out 
by  the  U-g  method.  Plots  of  flutter  speed  and  corresponding  reduced  frequency 
versus  alrfoll-alr  mass  ratio  are  shown  In  Figure  15.  The  values  assumed  for  the 
aeroelastlc  parameters  were:  position  of  elastic  axis  a^  * -0.5;  position  of 
mass  center  xa  * 0.25;  plunge  to  pitch  frequency  ratio  w^/w  “0.2;  radius  of 
gyration  r^  ■ 0.5;  and  damping  parameters  ^ " 0 anc*  !a  " 

The  results  obtained  from  the  relaxation  method  (UTRANS2)  and  the  Indicia! 
method  are, in  general, In  good  agreement.  The  curves  obtained  from  the  time 
integration  method,  however,  di  t'fer  from  the  other  two  sets  of  curves  to  some 
extent,  Such  discrepancy  may  be  due  to  the  differences  In  the  methods  and  the 
numerical  errors  In  computations.  It  may  be  noted  here  that  both  the  relaxation 
and  Indlclal  methods  are  based  on  time-linearized  equations  In  the  unsteady  aero- 
dynamic computations  and  the  provision  of  treating  the  shock  wave  motion  is  lacking. 
The  time  Integration  method  has  no  such  time-linearization  assumption  and  it  can 
efficiently  treat  shock  movement.  Thus  the  presence  of  shock  for  the  Mach  number 
considered  might  have  also  caused  the  discrepancy  In  these  flutter  results. 

In  the  response  analysis,  Equation  20  is  the  governing  aeroelastlc 
equation  of  the  system.  Reduced  frequency  was  assumed  as  kc  * 0.1.  To  start 
with,  the  aerodynamic  equation  was  integrated  in  time  for  six  cycles  by  forcing  a 
sinusoidal  pitching  motion  with  amplitude  of  0.01  radians.  In  the  sixth  cycle 
the  forced  responses  in  LTRAN2  became  almost  periodic.  From  this  point  onwards 
free  motion  was  considered  by  simultaneously  integrating  the  aerodynamic  and 
structural  equations.  Free  motion  was  started  with  initial  conditions  with 
f-(0)  - 0,  £'(0)  » 0,  a(0)  “ 0,  and  «'(0)  - 0.01. 
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Figure  15,  Effect  of  Airfoil-Air  Mass  Ratio  on  Flutter  Speed  for  NACA 
64A006  Airfoil  at  M - 0.85. 
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Stable,  neutrally  stable,  and  unstable  responses  were  obtained  by  parametrically 

varying  the  alrfoll-alr  mass  ratio  y . A neutrally  stable  response  was  obtained  by 

considering  a point  on  the  flutter  speed  curve  obtained  by  the  time  Integration 

method  In  Figure  15  at  k ■ 0.1.  The  corresponding  values  of  u and  U*  thus 

c 

obtained  were  equal  to  252.5  and  5.5,  respectl vely.  These  values  along  with  the 
other  aeroelastlc  parameters  (ah  « -.5,  x ■ 0.25,  u>K/u  ■ 0.2,  r ■ 0.5,  k » v.l) 
were  substituted  Into  Equation  20  and  the  response  analysis  was  carried  out. 

The  responses  for  plunging  displacement  and  corresponding  lifting  force  obtained 
from  LTRAN2  are  shown  In  Figure  16.  After  some  small  Initial  disturbances  due  to 
the  initial  conditions,  the  response  curves  show  Ideal  periodic  behavior. 

The  corresponding  stable  response  curves  for  pitching  displacement  a and 
pitching  moment  cm  are  shown  In  Figure  17.  These  curves  are  shown  only  from  the 
second  cycle  of  the  free  motion.  As  compared  to  Figure  16,  there  were  quite 
pronounced  initial  disturbances  due  to  Initial  conditions  in  the  first  cycle.  How- 
ever, the  responses  became  almost  periodic  after  the  third  cycle  of  free  motion. 

One  possible  reason  for  the  pronounced  initial  disturbances  In  the  first 
cycle  of  the  free  pitching  response  Is  that  the  pitching  axis  Is  located  at  the 
quarter  chord.  The  moment  values  at  the  quarter  chord  are  very  small  and  the 
plunging  motion  Initially  dominates  the  response.  Hence, the  pitching  responses 
require  more  cycles  to  become  periodic  when  compared  to  the  plunging  responses. 

It  may  be  noted  that  for  the  case  of  flat  plate  pitching  about,  the  mid-chord, 
the  order  of  Initial  disturbances  was  the  same  for  both  plunging  and  pitching 
responses  (See  Figures  11  and  12).  This  may  be  due  to  the  fact  that  the  moment 
value  about  the  mid-chord  is  larger  than  that  about  the  quarter  chord. 

The  curves  for  neutrally  stable  responses  presented  in  Figures  16  and  17  show 
that  the  neutrally  stable  conditions  obtained  by  the  present  response  method  are  In 
good  agreement  with  those  obtained  by  the  U-g  method  based  on  the  time  integration 
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Figure  16.  Effect  of  Airfoil -Air  Mass  Ratio  (Altitude)  on  Plunging  Oisplaceeent  (e)  and  Lifting  Force 
Coefficient  (c. ) for  HACA  64AG06  Airfoil  at  M = 0.85. 


Effect  of  Airfoil -Air  Mass  Ratio  (Altitude)  on  Pitching  Rotation 
Pitching  Honsent  (c)  for  NACA  64A006  Airfoil  at  H =0.85. 


method.  On  the  other  hand  the  flutter  points  predicted  by  the  Indlclal  and  re- 
laxation methods  will  result  In  a slightly  diverging  case. 

The  response  results  for  the  neutrally  stable  condition  presented  In  Figures 
16  and  17  also  show  the  validity  of  the  principle  of  superposition  of  air  loads 
which  was  used  in  obtaining  the  flutter  results  by  the  ll-g  method.  In  deriving  the 
equations  for  the  simultaneous  Integration  procedure  It  was  not  necessary  to  use 
such  assumption. 

In  Figure  16,  the  response  curves  for  plunging  displacement  t and  lifting 
force  c£  for  stable  and  unstable  conditions  are  also  shown.  These  were  obtained  by 
changing  the  airfoil -air  mass  ratio  u.  For  unstable  response,  the  v>  value  assumed 
was  lot  less  than  that  corresponding  to  the  neutrally  stable  condition.  This 
point  Is  In  the  unstable  zone  of  Figure  15.  On  the  other  hand,  for  stable  response 
the  value  of  u assumed  was  1 0X  higher  than  that  corresponding  to  the  neutrally 
stable  condition.  This  point  Is  In  the  stable  zone  of  Figure  15.  Stable  and 
unstable  response  curves  obtained  for  pitching  displacement  a and  pitching  moment 
are  also  shown  In  Figure  17.  It  Is  seen  that  the  stable  and  unstable  conditions 
obtained  by  the  U-g  method  produced  converging  and  diverging  responses,  respectively. 

It  was  also  of  Interest  to  study  the  effect  of  alrfoll-alr  mass  ratio  u 
(altitude)  on  pe^k  amplitudes  of  response  curves.  In  this  analysis,  the  peak 
amplitudes  corresponding  to  the  fourth  cycle  of  free  motion  were  considered, 

Figure  18  shows  plots  of  the  ratio  of  fourth  cycle  peak  amplitude  to  the 
amplitude  of  the  neutrally  stable  curves  versus  the  ratio  of  alrfoll-alr  mass 
ratio  u to  the  w value  corresponding  to  the  neutrally  stable  point.  These  plots 
Include  four  curves  for  plunging  displacement,  lifting  force,  pitching  displace- 
ment and  pitching  moment,  respectively. 

From  Figure  18  It  may  be  observed  that  all  the  four  curves  have  similar  trends. 
The  curve  corresponding  to  pitching  rotation  Is  the  steepest  among  the  four  curves. 
There  Is  a sudden  change  In  the  slopes  near  the  flutter  point  for  all  curves. 
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Figure  18.  Effect  of  Airfoil-Air  Mass  Ratio  on  Peak  Amplitudes 
of  Responses  for  NACA  64A006  Airfoil  at  M ■ 0.85. 
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SECTION  V 


CONCLUDING  REMARKS 

The  purpose  of  this  research  was  to  study  the  aeroelastlc  response  behavior 
of  single  (plunging  or  pitching)  and  two  degree  of  freedom  airfoil  systems  In 
two-dimensional  small-disturbance  unsteady  transonic  flows.  The  aerodynamic 
computer  code  LTRAN2  which  Is  based  on  a fully  Implicit  finite  difference 
scheme  was  employed  to  obtain  the  aerodynamic  forces.  Response  behavior  was 
studied  for  a flat  plate  at  a subsonic  Mach  number  of  M ■ 0.7  and  a NACA  64A006 
airfoil  at  transonic  Mach  numbers  0.85  and  0,88.  The  Influence  of  the  alrfoll- 
alr  mass  ratio  u on  response  behaviors  was  studied  In  detail. 

As  a result  of  this  study  the  following  concluding  remarks  may  be  made: 

(1)  The  response  results  obtained  for  single  degree  of  freedom  systems  of  pitching 
and  plunging  for  flat  plates  at  M « 0.7  show  that  LTRAN2  (linear)  results  compare 
well  with  the  results  based  on  quasi  steady-state  aerodynamic  theory. 

(2)  Pitching  response  results  obtained  for  a NACA  64A006  airfoil  at  M * 0.88 
show  that  the  present  procedure  of  numerically  Integrating  the  structural  and 
transonic  aerodynamic  equations  gives  results  which  agree  with  those  obtained 
previously  by  Ball haus  and  Goorjlan  (Reference  4). 

(3)  Results  obtained  for  a flat  plate  with  two  degrees  of  freedom  (pitching  and 
plunging)  at  M * 0.7  Illustrate  that  neutrally  stable  conditions  obtained  by  the 
U -g  method  based  on  aerodynamic  coefficients  computed  by  the  time  integration 
method  (LTRAN2)  check  with  the  neutrally  stable  conditions  obtained  by  the  time- 
response  method.  It  was  also  Illustrated  that  the  neutrally  stable  conditions 
obtained  by  the  U-g  method  based  on  aerodynamic  coefficients  obtained  by  the 
Kernel  Function  method  compare  well  with  those  obtained  by  the  time-response  method, 
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(4)  Response  results  obtained  for  a NACA  64A006  airfoil  with  two  degrees  of 
freedom  (pitching  and  plunging)  at  the  transonic  Mach  number  M - 0.85  show  that 
neutrally  stable  conditions  obtained  by  the  U-g  method  based  on  aerodynamic 
coefficients  computed  by  the  time  Integration  method  (LTRAN2)  agree  well  with 
those  obtained  by  the  time  Integration  response  analysis.  However,  the  neutrally 
stable  conditions  obtained  by  the  U-g  method  based  on  the  relaxation  and  indlclal 
aerodynamic  coefficients  do  not  agree  so  well  with  that  obtained  by  the  time 
Integration  response  analysis,  This  lack  of  good  agreement  may  be  due  to  the 
assumption  of  time  linearization  In  the  Indlclal  and  relaxation  methods.  Further- 
more, the  provision  of  treating  the  shock  wave  motion  Is  lacking  In  the  relaxation 
and  the  Indlclal  programs. 

(5)  Good  comparison  between  the  neutrally  stable  conditions  given  by  the  U-g 
method  (based  on  time  Integration  aerodynamic  coefficients)  and  the  time- 
response  analysis  Indicate  that  the  principle  of  linear  superposition  of  airloads 
used  In  the  U-g  method  Is  valid  for  the  cases  analyzed. 

(6)  In  general  It  may  be  concluded  that  the  present  method  Is  accurate  for 
predicting  the  neutrally  stable  conditions  (flutter)  for  a two-degree-of-freedom 
system.  This  method  also  takes  Into  consideration  the  movement  of  the  shock 
whereas  It  Is  not  possible  to  do  so  In  the  U-g  method, 


(7)  Because  of  the  limitations  of  the  method  employed  In  the  computer  code 
LTRAN2 , only  low  reduced  frequencies  (k.  < 0.2)  were  considered  In  the  analysis. 

V 

Modifications  to  the  LTRAN2  code  are  required  for  consideration  of  higher 
reduced  frequencies. 

(8)  In  this  study  It  was  also  observed  that  M ■ 0.88  appears  to  be  practically 

the  highest  Mach  number  that  can  be  used  In  LTRAN2  for  a NACA  64A006.  Modifications 
In  the  basic  code  are  required  to  account  for  higher  Mach  numbers. 
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(9)  This  analysis  may  provide  a systematic  procedure  for  the  transonic  aero- 
elastic  response  analysis  of  airfoils.  The  present  study  may  also  be  used  as 
a comparative  basis  by  the  experimentalists  and  other  analysts  who  work  in  the 
same  field. 

(10)  Similar  studies  on  other  airfoils,  especially  supercritical  ones,  are 
needed. 

(11)  In  order  to  study  more  practical  cases  of  full-wing  aeroelastlc  responses, 
a transonic  computer  code  to  account  for  three  dimensional  unsteady  flow  Is 
required. 
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